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ABSTRACT 

We  have  imaged  a  ~1  deg2  field  centered  on  the  known  Galactic  supernova  remnant  (SNR)  G1 1.2— 0.3  at  74, 
330,  and  1465  MHz  with  the  Very  Large  Array  radio  telescope  and  235  MHz  with  the  Giant  Metrewave  Radio 
Telescope.  The  235,  330,  and  1465  MHz  data  have  a  resolution  of  25",  while  the  74  MHz  data  have  a  resolution  of 
~100".  The  addition  of  this  low-frequency  data  has  allowed  us  to  confirm  the  previously  reported  low-frequency 
turnover  in  the  radio  continuum  spectra  of  the  two  known  SNRs  in  the  field,  Gil. 2— 0.3  and  G11.4— 0.1,  with 
unprecedented  precision.  Such  low-frequency  turnovers  are  believed  to  arise  from  free-free  absorption  in  ionized 
thermal  gas  along  the  lines  of  sight  to  the  SNRs.  Our  data  suggest  that  the  74  MHz  optical  depths  of  the  absorbing 
gas  is  0.56  and  1.1  for  G1 1.2— 0.3  and  Gil. 4— 0.1,  respectively.  In  addition  to  adding  much  needed  low-frequency 
integrated  flux  measurements  for  two  known  SNRs,  we  have  also  detected  three  new  SNRs:  Gil. 15— 0.71, 
Gil. 03— 0.05,  and  Gil. 18+0. 11.  These  new  SNRs  have  integrated  spectral  indices  between  —0.44  and  —0.80. 
Because  of  confusion  with  thermal  sources,  the  high  resolution  (compared  with  previous  Galactic  radio  frequency 
surveys)  and  surface  brightness  sensitivity  of  our  observations  have  been  essential  to  the  identification  of  these  new 
SNRs.  With  this  study  we  have  more  than  doubled  the  number  of  SNRs  within  just  a  ~1  deg2  field  of  view  in  the 
inner  Galactic  plane.  This  result  suggests  that  future  low-frequency  observations  of  the  Galactic  plane  of  similar 
quality  may  go  a  long  way  toward  alleviating  the  long-recognized  incompleteness  of  Galactic  SNR  catalogs. 
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1.  INTRODUCTION 

Supernova  explosions  have  a  profound  effect  on  the  mor¬ 
phology,  kinematics,  and  ionization  balance  of  galaxies,  and 
possibly  trigger  new  generations  of  star  formation.  However, 
based  on  statistical  studies  of  the  incompleteness  of  Galactic 
supernova  remnant  (SNR)  catalogs  and  predictions  of  the 
Galactic  SNR  rate,  there  should  be  many  more  SNRs  in  our 
Galaxy  (~  1 000;  Li  et  al.  1991;  Tammann,  Loeffler,  & 
Schroeder  1994;  Case  &  Bhattacharya  1998)  than  are  cur¬ 
rently  known  (~230:  Green  2002).  This  paucity  is  likely  due 
in  part  to  selection  effects  acting  against  the  discovery  of  the 
more  mature,  faint,  extended  remnants,  as  well  as  of  the  very 
young,  small  remnants  as  a  result  of  poor  sensitivity  and 
spatial  resolution  at  low  frequencies  where  SNRs  are  brightest 
(see  Green  1991).  These  missing  remnants  are  thought  to  be 
concentrated  toward  the  inner  Galaxy,  where  the  diffuse 
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nonthermal  Galactic  plane  emission,  coupled  with  the  thermal 
emission  from  H  n  regions,  causes  the  most  confusion.  A 
complete  census  of  Galactic  SNRs  is  essential  to  under¬ 
standing  the  star  formation  history  of  our  Galaxy. 

Recent  advances  in  both  the  instrumentation  and  software 
for  low  radio  frequency  observations  make  sensitive,  high- 
resolution  imaging  in  this  relatively  unexplored  wave  band 
possible  for  the  first  time  (see,  e.g.,  LaRosa  et  al.  2000).  As  a 
test  case  for  a  larger  fully  sampled  mosaic  of  the  inner 
Galactic  plane  from  l  =  +4°  to  +20°,  we  imaged  the  region 
around  the  known  SNR  G1 1.2— 0.3  with  the  National  Radio 
Astronomy  Observatory  (NRAO)9  Very  Large  Array  (VLA)  at 
74  and  330  MHz  using  data  from  multiple  configurations.  We 
have  also  obtained  a  VLA  mosaic  at  1465  MHz  of  this  region 
along  with  235  MHz  data  from  the  Giant  Metrewave  Radio 
Telescope  (GMRT).  The  final  region  of  overlap  of  these 
different  data  sets  covers  a  ~1  deg2  field  of  view  in  which  we 
detect  two  known  SNRs,  G1 1.2— 0.3  and  Gil. 4— 0.1,  and 
identify  three  new  SNRs.  The  results  from  these  observations 
are  presented  below. 

2.  OBSERVATIONS 

We  have  carried  out  observations  of  the  Galactic  plane 
centered  on  the  known  SNR  Gil. 2— 0.3  at  four  frequencies: 
1465,  330,  235,  and  74  MHz.  The  1465,  330,  and  74  MHz  data 
were  obtained  with  the  VLA  in  multiple  configurations,  while 
the  235  MHz  data  were  obtained  from  the  GMRT.  The  FWHP 
primary  beams  of  the  single  pointing  330,  235,  and  74  MHz 
data  are  2°.  5,  1  ? 5,  and  11(5,  respectively.  Data  at  all  four 
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frequencies  were  acquired  in  spectral  line  mode  to  aid  in  radio 
frequency  interference  (RFI)  excision  and  minimize  bandwidth 
smearing.  The  parameters  for  these  observations  are  listed  in 
Table  1 .  The  data  from  all  four  frequencies  were  reduced  and 
imaged  using  the  wide-field  imaging  capabilities  within  the 
NRAO  Astronomical  Image  Processing  Software  (AIPS) 
package.  More  information  on  wide-field  imaging  techniques 
can  be  found  in  Cornwell  &  Perley  (1992).  Additional  details  of 
wide-field  Galactic  plane  imaging  can  be  found  in  LaRosa  et  al. 
(2000). 

The  initial  antenna-based  gain  and  phase  corrections  for  the 
235,  330,  and  1465  MHz  data  were  calculated  in  the  usual 
way  from  observations  of  the  phase  calibrator  J1833— 210  and 
either  3C286  or  3C48  for  flux  calibration.  Unfortunately,  there 
are  no  pointlike  calibrators  at  74  MHz  that  are  strong  enough 
to  dominate  the  total  flux  in  the  large  74  MHz  primary  beam 
(11?  5).  For  this  reason,  the  flux  scale  and  a  single  instrumental 
phase  at  74  MHz  were  estimated  from  observations  of  Cygnus 
A  using  a  publicly  available  model10  of  its  source  structure. 
Since  Cygnus  A  is  far  from  the  target  field,  initial  estimates  for 
the  phase  as  a  function  of  time,  as  well  as  positional 
information  for  the  target  field,  were  obtained  by  using  the 
VLA  330  MHz  image  as  the  model  in  the  first  round  of 
74  MHz  phase-only  self-calibration  (see  Kassim  et  al.  2004, 
for  further  details  specific  to  74  MHz  VLA  data  reduction). 
At  all  four  frequencies,  following  the  initial  calibrator-based 
corrections,  several  iterations  of  self-calibration  were  applied 
to  data  from  each  configuration  or  observation  listed  in  Table  1 
to  improve  the  dynamic  range.  The  calibrated  visibility  data 
from  all  of  the  available  observations  for  each  frequency  were 
then  concatenated  into  a  single  U—  V  data  set,  after  which  one 
further  amplitude  and  phase  calibration  was  performed. 

The  1465  MHz  VLA  image  is  composed  of  a  three  pointing 
mosaic,  with  the  pointings  arranged  in  a  triangle  with 
Gil. 2— 0.3  near  its  center.  The  FWHP  primary  beam  of  each 
pointing  is  0 .'  5 .  The  three  pointings  were  calibrated  and 


10  See  http://lofar.nrl.navy.mil/pubs/tutoriaWLAmodels/models.html. 


imaged  separately  and  then  combined  with  the  linear 
mosaicking  task  FLATN  within  AIPS.  During  the  FLATN 
process  the  individual  pointings  were  corrected  for  the  re¬ 
sponse  of  the  primary  beam  and  then  masked  at  the  25%  power 
point.  The  final  1465  MHz  image  covers  a  ~1  deg2  area.  After 
correction  for  the  primary  beam  response,  the  330,  235,  and  74 
MHz  images  were  subimaged  to  the  same  field  of  view. 

To  correct  for  the  gain  compression  that  the  GMRT  suffers 
when  observing  toward  bright  sources  (i.e.,  the  inner  Galactic 
plane),  an  additional  gain  correction  factor  must  be  applied  to 
the  235  MHz  data  after  imaging.  This  factor  can  be  estimated 
at  the  time  of  the  observations  by  comparing  the  total  power 
data  for  the  calibrators  with  and  without  the  automatic  level 
attenuators  on  for  each  antenna.  At  the  time  of  the  second 
GMRT  235  MHz  observation  (see  Table  1),  this  empirical 
factor  had  an  average  value  of  4  (this  factor  is  not  available  for 
the  first  GMRT  observation).  However,  the  integrated  flux 
density  measurements  of  the  known  SNRs  (Gil. 2— 0.3  and 
G1 1.4— 0.1),  after  multiplying  the  final  image  by  this  cor¬ 
rection  factor,  were  excessively  high  given  their  known  spec¬ 
tral  indices.  This  discrepancy  may  be  due  in  part  to  the 
amplitude  self-calibrations  that  were  applied  to  the  235  MHz 
data  before  the  gain  correction  factor,  or  it  is  possible  that 
the  correction  factor  for  the  first  observing  run  may  be  differ¬ 
ent  from  that  of  the  second.  In  order  to  better  calibrate  the 
235  MHz  gain,  we  have  used  the  spectral  index  between  330 
and  1465  MHz  at  each  pixel  for  Gil. 2— 0.3  and  Gil. 4— 0.1  to 
create  an  image  of  the  expected  flux  density  at  235  MHz  for 
each  pixel.  This  approach  is  justified  since  the  continuum 
spectra  of  SNRs  above  150  MHz  are  well  characterized  by 
a  constant  power-law  spectrum  (see,  e.g.,  Kassim  1989a). 
Comparison  of  the  expected  235  MHz  image  with  the 
observed  (uncorrected)  image  yielded  an  average  correction 
factor  of  2.7.  The  final  235  MHz  image  was  multiplied  by  2.7 
before  primary  beam  correction.  The  corrected  235  MHz 
image  was  then  convolved  to  a  resolution  of  25". 

For  the  1465  and  330  MHz  VLA  data,  care  was  taken  to 
match  the  U—  V  coverage  so  that  each  composite  data  set  had 
an  intrinsic  resolution  of  approximately  25".  The  images  were 


TABLE  1 

Gil. 2  Field  Observational  Parameters 


Bandwidth  Timea 

Date  Instrument  Config.  (MHz)  (hr) 


74/330  MHz  Parameters 

2001  Jan  13 .  VLA  A  . .  73  4.4 

2001  Mar  1 .  VLA  B  1.5/3  1.7 

2002  Aug  20 .  VLA  B  ..712  1.5 

2001  Aug  28 .  VLA  C  1.5/3  0.75 

2001  Sep  28 .  VLA  D  1.5/3  1.7 


235  MHz  Parameters 

2001  Sep  23b .  GRMT  ...  16  1.5 

2002  May  05 .  GMRT  ...  8  0.5 


1465  MHz  Parameters 

2001  Jun  27 .  VLA  CnB  25  0.3 

2001  Aug  03 .  VLA  C  25  0.3 

2001  Sep  24 .  VLA  DnC  25  0.3 


a  Approximate  final  time  on  source. 
b  Data  contributed  by  P.  Rao. 
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then  convolved  to  exactly  25",  but  the  convolution  was 
minimal.  The  1465  and  330  MHz  data  are  not  sensitive  to 
structures  larger  than  ~15'.  The  VLA  74  MHz  data  have  a 
resolution  of  145"  x  92".  The  74  MHz  image  does  not  include 
A-configuration  data  because  exceptionally  poor  ionospheric 
conditions  rendered  them  uncalibratable.  Since  the  74  MHz 
resolution  only  marginally  resolves  the  SNRs,  and,  as 
described  above,  the  amplitude  calibration  of  the  235  MHz 
data  is  a  bit  uncertain,  data  from  these  two  frequencies  were 
only  used  for  integrated  flux  measurements. 

3.  RESULTS 

Figures  la  and  lb  show  our  VLA  1465  and  330  MHz 
images  with  25"  resolution  of  a  ~1  deg2  area  centered  at 
G1 1.2— 0.3.  The  rms  noise  in  the  1465  MHz  mosaicked  image 
increases  as  a  function  of  distance  from  the  center  of  the 
image,  where  it  is  ~1  mJy  beam-1,  to  about  2  mJy  beam-1 
near  the  edge.  The  rms  noise  level  in  the  330  MHz  image  is 
~3  mJy  beam-1,  while  the  (peak)  dynamic  range  in  these  two 
images  are  470  and  340,  respectively.  The  known  SNRs 
Gil. 2— 0.3  and  G11.4— 0.1  are  labeled  on  Figure  la  for 
reference,  along  with  the  locations  of  the  three  new  SNR 
candidates,  designated  A,  B,  and  C.  Four  criteria  were  used  to 
identify  the  new  SNR  candidates:  (1)  an  extended  region  of 
emission  needed  to  be  present  at  both  1465  and  330  MHz; 
(2)  the  extended  emission  needed  to  have  a  nonthermal  spectral 
index  typical  of  SNRs  (about  —0.1  to  —0.9);  (3)  the  extended 
emission  needed  to  have  a  shell-like  (or  partial  shell-like) 
morphology;  and  (4)  the  positions  of  sources  meeting  the  first 
three  criteria  were  then  searched  with  both  NED  and  S1MBAD 
to  be  sure  that  they  are  not  associated  with  any  previously 
identified  sources  of  nonthermal  emission.  The  Galactic 
coordinates  for  the  new  SNR  candidates  meeting  these  criteria 
are  A:  Gil. 15-0.71,  B:  Gll.03-0.05,  and  C:  Gil. 18+0.11. 
These  positions  have  been  taken  from  the  approximate  center 
of  the  SNR  shells,  and  for  the  purpose  of  this  paper  two 
significant  digits  (instead  of  the  usual  one)  are  used  to  help 
distinguish  the  new  candidate  SNRs  from  the  two  known 
remnants.  Given  the  strong  morphological  and  spectral 
evidence  that  is  presented  in  more  detail  below,  we  will 
henceforth  refer  to  these  sources  as  “new  SNRs,”  rather  than 
“SNR  candidates.” 

Integrated  flux  density  measurements  were  made  for  the  four 
frequencies  in  our  study  for  all  five  SNRs.  Since  the  images 
have  different  noise  levels  and  several  of  the  sources  have 
regions  of  “breakout”  morphology  that  tend  to  gradually  fade 
into  the  noise  (i.e.,  not  sharp  shells),  it  was  difficult  to  determine 
where  the  boundaries  of  the  SNRs  should  be  drawn.  Therefore, 
in  order  to  get  comparable  flux  density  estimates,  the  4  o 
(12  mJy  beam-1)  contour  of  the  330  MHz  image  (see  Fig.  lb) 
was  used  to  define  the  SNRs  for  all  of  the  integrated  flux 
measurements.  The  errors  on  these  flux  density  measurements 
were  estimated  by  (number  of  independent  beams)1/2  x  3<rat 
each  frequency.  Upper  limits  for  nondetections  were  estimated 
by  (number  of  independent  beams)1/2  x  5  a.  The  rms  noise  of 
the  235  and  74  MHz  data  are  8  and  150  mJy  beam-1,  respec¬ 
tively.  When  necessary,  an  extra  contribution  from  the 
background  was  subtracted  from  the  integrated  flux  density. 
This  contribution  was  estimated  by  integrating  the  flux  in  a 
nearby  source-free  region  with  the  same  approximate  area  as 
the  source.  The  background  level  correction  was  only 
necessary  for  the  74  and  235  MHz  images.  In  addition,  the 
error  estimates  for  the  235  MHz  flux  densities  have  had  an 
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Fig.  1. — (a)  VLA  1465  MHz  three-pointing  mosaic  centered  on  the  known 
SNR  Gil. 2— 0.3  with  25"  resolution.  The  contour  levels  are  4,  10,  20,  50, 
100,  200,  and  400  mJy  beam-1.  The  rms  noise  of  this  image  is  —  I  mJy 
beam  1  but  varies  somewhat  as  a  function  of  position  as  a  result  of  the 
mosaicking  procedure.  The  sources  identified  with  the  dotted  boxes  and 
labeled  A,  B,  and  C  are  the  three  new  SNR  candidates,  (b)  VLA  image  at 
330  MHz  centered  on  SNR  G1 1 .2—0.3  with  25"  resolution.  The  contour  levels 
are  10,  20,  50,  100,  200,  400,  and  700  mJy  beam  *.  The  rms  noise  of  this 
image  is  —3  mJy  beam-1.  The  region  outside  of  the  1465  MHz  mosaic  has 
also  been  blanked  on  the  330  MHz  image  for  ease  of  comparison. 

additional  5%  error  added  to  account  for  the  extra  uncertainty  in 
the  235  MHz  gain  calibration  described  above. 

Flux  densities  from  the  literature  have  also  been  included  in 
our  analysis  for  all  cases  in  which  an  error  estimate  is 
available  and  the  error  is  less  than  20%.  Figures  2  and  3  show 
the  integrated  radio  continuum  spectra  for  the  two  known 
SNRs  in  our  field,  Gil. 2— 0.3  and  Gil. 4— 0.1.  Figures  4-6 
show  the  integrated  radio  continuum  spectra  for  the  new  SNRs 
A  (Gil. 15-0.71),  B  (Gll.03-0.05),  and  C  (Gil. 18+0.11). 
Further  details  for  each  SNR  are  described  in  §  3.1. 
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1  2  3  4  5 

log  (Frequency  [MHz]) 

Fig.  2. — Radio  continuum  spectrum  for  SNR  Gil. 2— 0.3.  The  best-fit 
spectral  index  using  eq.  (1)  is  —0.570  ±  0.002.  The  hexagon  symbols  are  from 
the  current  work,  the  triangles  are  from  Kassim  (1989b,  1992),  the  single 
square  is  from  Kovalenko,  Pynzar,  &  Udaltsov  (1994),  and  the  crosses  are 
from  Kothes  &  Reich  (2001).  To  obtain  the  spectral  index  for  only  the  shell  of 
the  SNR,  1  Jy  has  been  subtracted  from  the  integrated  flux  measurements  for 
data  with  v  >  200  MHz,  assuming  that  ap  ~  0.0. 

For  each  SNR  the  integrated  spectral  index  was  calculated 
using  a  weighted  least-squares  fit  using  the  data  points  shown 
in  Figures  2-6  (except  for  the  74  MHz  upper  limits).  Because 
of  the  method  described  above  to  calibrate  the  235  MHz  gain, 
the  data  points  at  235  MHz  were  not  used  to  fit  the  spectral 
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Fig.  4. — Radio  continuum  spectrum  for  new  SNR  A:  G1 1.15—0.71.  The 
best-fit  power-law  spectral  index  is  a  =  —0.82  ±0.01.  The  hexagon  symbols 
are  from  the  current  work,  while  the  integrated  flux  data  point  marked  by  the 
cross  is  from  Reich  et  al.  (2001).  The  flux  density  marked  by  the  open 
hexagon  at  74  MHz  is  an  upper  limit  and  was  not  used  in  the  fit. 

index  for  Gil. 2— 0.3  and  G1 1.4— 0.1  (i.e.,  since  these  235  MHz 
flux  densities  were  used  to  set  the  235  MHz  gain  scale).  The 
data  for  G1 1.2— 0.3  and  G1 1.4— 0.1  were  fitted  to  the  equation 

Sv  =  51000  (lOOOMHz)  6XP  [  ■Tl00°  (lOOOMHz)  J  ’  ( 1 } 
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Fig.  3. — Radio  continuum  spectrum  for  SNR  G1 1.4— 0.1.  The  best-fit 
spectral  index  using  eq.  (1)  is  —0.59  ±0.01.  The  hexagons  are  from  the 
current  work,  the  triangles  are  from  Kassim  (1989b)  and  Kassim  (1992),  and 
the  crosses  are  from  Reich,  Reich,  &  Furst  (1997)  and  Reich  et  al.  (2001). 
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Fig.  5. — Radio  continuum  spectrum  for  new  SNR  B:  G1 1.03—0.05.  The 
best-fit  power-law  spectral  index  is  —0.60  ±  0.01.  The  hexagons  are  from  the 
current  work.  The  flux  density  marked  by  the  open  hexagon  at  74  MHz  is  an 
upper  limit  and  was  not  used  in  the  fit. 
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Fig.  6. — Radio  continuum  spectrum  for  new  SNR  C:  Gil. 18+0.11.  The 
best-fit  power-law  spectral  index  is  —0.44  ±  0.01.  The  hexagons  are  from  the 
current  work.  The  flux  density  marked  by  the  open  hexagon  at  74  MHz  is  an 
upper  limit  and  was  not  used  in  the  fit. 

where  S[000  and  Tiooo  are  the  flux  density  and  optical  depth  at 
a  fiducial  frequency  of  1000  MHz,  respectively,  and  a  is  the 
integrated  spectral  index  (Dulk  &  Slee  1975;  Kassim  1989a). 
This  equation  assumes  a  standard  nonthermal  constant  power- 
law  spectrum  and  allows  for  a  thermal  absorption  turnover  at 
lower  frequencies  (see,  e.g.,  Dulk  &  Slee  1975;  Kassim  1989a). 
The  free-free  continuum  optical  depth  at  other  frequencies  can 
be  estimated  from  t„  =  riooo(^/1000  MHz)-2'1.  The  new 
SNRs  were  deemed  to  have  too  few  data  points  to  adequately 
constrain  equation  (1);  instead  they  were  fitted  with  a  pure 
power  law  (i.e.,S„  a  v+a).  The  solid  lines  in  Figures  2-6  show 
the  best  fits  for  each  SNR. 

Figures  la-le  show  spatially  resolved  spectral  index  maps 
toward  the  five  SNRs  between  330  and  1465  MHz.  Before 
calculating  the  spectral  index,  the  two  images  were  masked  at 
the  4  a  level  (12  and  4  mJy  beam-1,  respectively).  Because  of 
the  low-level  calibration  and  positional  errors  that  are  inherent 
in  these  low-frequency  images,  small-scale  spectral  index  vari¬ 
ations  less  than  ~  0 . 1  should  be  viewed  skeptically.  However, 
these  maps  are  quite  useful  for  identifying  regions  of  thermal 
(positive  a)  versus  nonthermal  emission  (negative  a). 

3.1.  Individual  Sources 
3.1.1.  SNR  Gil. 2-0. 3 

The  SNR  Gil. 2— 0.3  is  thought  to  be  coincident  with  the 
historical  supernova  of  A.D.  386  (Clark  &  Stephenson  1977). 
Although  there  is  still  some  ongoing  debate  over  the  reality  of 
this  association,  it  is  clear  that  G1 1.2— 0.3  is  quite  young 
(Roberts  et  al.  2003).  The  distance  to  this  remnant  has  been 
estimated  from  H  i  absorption  measurements  to  be  ~5  kpc 
(Radhakrishnan  et  al.  1972;  Green  et  al.  1988).  In  the  radio 
Gil. 2— 0.3  has  a  steep  spectrum  radio  shell  with  a  signifi¬ 
cantly  flatter  spectrum  central  core  (see  Fig.  la;  Tam,  Roberts, 
&  Kaspi  2002;  Morsi  &  Reich  1987;  Kothes  &  Reich  2001). 
The  soft  X-ray  morphology  of  Gil. 2— 0.3  follows  that  of  the 


radio  shell,  while  there  is  a  nonthermal,  hard  X-ray  core 
consistent  with  a  pulsar  wind  nebula  (PWN;  Vasisht  et  al. 
1996).  Subsequent  X-ray  observations  have  confirmed  the 
presence  of  a  65  ms  X-ray  pulsar  near  the  center  of  the  remnant 
(Torii  et  al.  1997;  Kaspi  et  al.  2001;  Roberts  et  al.  2003). 

Kothes  &  Reich  (2001)  have  studied  the  morphology  of 
G11.2— 0.3  between  10  and  32  GHz  using  the  Effelsberg  radio 
telescope,  and  they  find  that  the  data  are  consistent  with  a  radio 
shell  spectral  index  of  as  ~  —0.57  and  a  PWN  spectral  index 
of  ap  ~  0.0,  with  the  PWN  contributing  about  1  Jy  of  flux  to 
the  integrated  flux  density  of  the  remnant  (~6%  at  1 .4  GHz). 
A  higher  resolution  study  of  the  spectral  index  variations  of 
Gil. 2— 0.3  between  1465  and  4865  MHz  using  archival  VLA 
data  has  been  carried  out  by  Tam  et  al.  (2002)  to  study  the 
interaction  of  the  PWN  with  the  shell.  These  authors  find  that 
—0.56  and  that  the  PWN  region  has  ap~  —0.25  and 
contributes  about  1  Jy  of  flux  density  at  1465  MHz.  However, 
this  analysis  may  suffer  from  missing  short  spacings,  along 
with  confusion  from  the  SNR  shell  emission. 

To  account  for  the  flatter  spectrum  of  the  core,  1  Jy  has  been 
subtracted  from  the  integrated  flux  measurements  shown  in 
Figure  2  for  data  with  v  >  200  MHz,  assuming  that  a  P  ~  0.0. 
We  have  not  subtracted  1  Jy  below  200  MHz,  under  the 
assumption  that  the  PWN  could  suffer  from  self-absorption  at 
some  frequency,  as  well  as  the  extrinsic  free-free  absorption 
observed  for  the  SNR  as  a  whole.  However,  without  more 
information  about  the  spectral  properties  of  the  PWN,  this 
choice  of  cutoff  frequency  is  somewhat  arbitrary.  It  is  clear 
from  the  fit  shown  in  Figure  2  that  there  is  still  a  small  excess 
of  flux  at  higher  frequencies,  suggesting  that  ap  is  actually 
somewhat  inverted  (i.e.,  ap  >  0.0).  A  more  detailed  study  of 
the  spectral  characteristics  of  the  PWN  will  be  presented  in 
Tam  et  al.  (2004),  based  on  the  new  VLA  330  MHz  image  and 
1465  MHz  mosaicked  image  presented  here,  along  with  new 
4835  and  8600  MHz  mosaics  of  this  region. 

The  integrated  continuum  spectrum  shown  in  Figure  2  for 
Gil. 2— 0.3  shows  a  sharp  downturn  for  v  <  200  MHz.  This 
spectral  behavior  was  also  noted  by  Kassim  (1989a)  in  a 
survey  to  determine  the  spectral  properties  of  SNR  at  low 
frequencies,  and  it  will  be  discussed  further  in  §  4.2.  The  best- 
fit  parameters  to  the  data  shown  in  Figure  2  (after  subtracting 
1  Jy  as  described  above)  are  a  =  —0.570  ±  0.002  for  the  shell 
and  an  optical  depth  at  74  MHz  of  r74  ~  0.56.  The  spatially 
resolved  spectral  index  of  Gil. 2— 0.3  is  shown  in  Figure  la. 
The  flatter  spectrum  of  the  PWN  is  clearly  visible  in  this 
image. 

The  nature  of  the  weak,  extended  (~15')  linear  feature  that 
lies  tangent  to  the  southern  edge  of  Gil. 2— 0.3  (see  Figs,  la, 
lb,  and  la)  is  unknown.  The  radio  spectrum  of  this  extended 
emission  is  fairly  steep,  with  spectral  indices  of  about  —0.5. 
Moreover,  there  are  no  extended  infrared  sources  near  this 
linear  structure  (see  §  4.1).  Together,  these  facts  suggest  that 
the  emission  is  nonthermal  in  origin.  Although  this  source 
meets  three  of  the  four  criteria  described  in  §  3  to  be  an  SNR 
candidate,  it  does  not  appear  to  have  shell-like  morphology. 
Instead,  it  appears  to  be  both  spectrally  and  morphologically 
similar  to  the  linear  nonthermal  radio  filaments  (NTFs) 
observed  toward  the  Galactic  center  region  (see,  for  example, 
LaRosa  et  al.  2000).  However,  to  our  knowledge,  NTFs  have 
not  previously  been  observed  outside  the  immediate  Galactic 
center  region  (i.e.,  within  ~1°).  The  nature  of  the  Galactic 
center  NTFs  is  also  unclear,  although  many  of  them  are  ap¬ 
proximately  perpendicular  to  the  Galactic  plane,  suggesting 
that  they  may  be  associated  with  magnetic  field  enhancements 
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Fig.  7. — Spectral  index  maps  between  330  and  1465  MHz  for  ( a )  Gil. 2— 0.3,  ( b )  G1 1.4— 0.1,  (c)  new  SNR  A:  Gil.  15— 0.71,  ( d )  new  SNR  B:  G1 1.03—0.05,  and 
( e )  new  SNR  C:  G1 1. 18+0.1 1.  The  1465  MHz  contours  displayed  in  Fig.  1  a  are  also  shown.  The  330  and  1465  MHz  continuum  images  were  masked  at  12  and  4  mJy 
beam-1,  respectively,  before  the  spectral  indices  were  calculated.  Note  that  the  color  scale  changes  from  image  to  image. 


(see,  e.g.,  Morris  &  Serabyn  1996).  Higher  resolution  and 
sensitivity  observations  of  the  G1 1.2— 0.3  linear  feature  may 
help  to  shed  light  on  this  interesting  source. 

3.1.2.  SNR  Gil. 4-0.1 

G1 1.4— 0.1  was  first  identified  as  an  SNR  by  Caswell  et  al. 
(1975).  The  distance  to  this  SNR  is  currently  unknown.  The 
best-fit  parameters  for  the  integrated  continuum  spectrum 
shown  in  Figure  3  are  a  =  —0.59  ±  0.01  and  an  optical  depth 
at  74  MHz  of  T74  ~  1 . 1 .  The  spatially  resolved  spectral  index 
map  for  G1 1.4— 0.1  shown  in  Figure  lb  suggests  the  presence 
of  a  number  of  flatter  (a~  —0.4)  than  average  (a~—  0.6) 
spectral  index  filaments  running  through  the  SNR.  Similar 
spectral  flattening  toward  filamentary  structures  has  recently 
been  observed  toward  a  number  of  mature  SNRs,  including 
W28  and  W49B  (Dubner  et  al.  2000;  Moffett  &  Reynolds 
1994).  Such  flattening  may  be  indicative  of  first-order  Fermi 
shock  acceleration  at  the  sites  of  strong  shocks  (see,  e.g.,  Jun 
&  Jones  1999). 

Based  on  63"  x  44"  resolution  1465  MHz  VLA  images, 
Dubner  et  al.  (1993)  suggest  that  Gil. 4— 0.1  could  be  com¬ 
posed  of  two  superposed  remnants.  From  the  330/1465  MHz 
spectral  index  map  shown  in  Figure  lb,  it  seems  likely  that 


G1 1.4— 0.1  is  composed  of  a  single  nonthermal  structure.  The 
point  source  near  R.A.  =  18h10m36s,  deck  =  -19°08'00" 
(J2000.0)  in  Figures  la,  lb,  and  lb  has  a  flat  radio  spectrum 
( a  ~  +0.1),  and,  since  it  does  not  have  an  infrared  counterpart, 
it  is  unlikely  to  be  thermal  in  origin  (see  §  4.1)  and  is  most 
likely  an  unrelated  extragalactic  source.  It  could  also  be  a 
pulsar,  but  this  region  of  the  plane  was  observed  in  the  Parkes 
Multibeam  Pulsar  Survey  (Morris  et  al.  2002)  and  it  seems 
unlikely  that  a  pulsar  with  such  a  strong  radio  continuum 
would  have  gone  undetected.  While  strong  to  be  a  heretofore 
undetected  pulsar,  the  point  source  only  contributes  a  few  tens 
of  millijanskys  to  the  integrated  SNR  flux  density,  so  no  effort 
was  made  to  subtract  it  from  the  integrated  flux  measurements. 

3.1.3.  New  SNR  A:  Gil. 15-0.71 

The  integrated  continuum  spectrum  for  the  new  SNR  A  lo¬ 
cated  at  Gil. 15—  0.71  is  shown  in  Figure  4,  while  a  map  of  its 
spatially  resolved  spectral  index  between  330  and  1465  MHz 
is  shown  in  Figure  7c.  At  330  and  1465  MHz  Gil.  15— 0.71 
appears  as  a  rim-brightened  partial  shell  on  its  northern  end. 
Gil. 15— 0.71  has  also  been  weakly  detected  in  the  2695  MHz 
Bonn  survey  (Reich  et  al.  2001;  this  data  point  is  included 
on  Fig.  4)  but  has  not  previously  been  spatially  resolved  or 
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recognized  as  a  nonthermal  source.  Its  integrated  spectral 
index  is  quite  steep,  with  a  =  —0.82  ±  0.01.  The  boundary  of 
this  new  SNR  is  particularly  difficult  to  determine  given  its  low 
surface  brightness  and  the  increased  noise  toward  the  edge  of 
the  1465  MHz  image  (it  is  unclear  whether  the  emission  to¬ 
ward  the  southwestern  edge  of  region  A  in  Fig.  la  is  real). 
However,  the  goodness  of  the  spectral  index  fit  shown  in 
Figure  4,  which  includes  the  single-dish  2695  MHz  data  point, 
suggests  that  we  are  not  missing  very  much  of  the  remnant’s 
total  flux.  Nevertheless,  higher  sensitivity  pointed  observa¬ 
tions  of  this  source  are  needed  to  determine  its  full  extent  and 
morphology. 

3.1.4.  New  SNR  B:  Gll.03-0.05 

New  SNR  B  at  G1 1.03—0.05  appears  as  a  complete  shell  at 
330  MHz  and  a  partial  shell  at  1465  MHz  (see  Figs,  la  and  lb). 
A  plot  of  its  integrated  continuum  spectrum  is  shown  in 
Figure  5,  while  the  morphology  of  its  spatially  resolved 
330/1465  MHz  spectral  index  is  shown  in  Figure  Id. 
G1 1.03—0.05  has  an  integrated  spectral  index  of  — 0.60  ± 
0.01.  The  double  source  at  R.A.  =  18h09m43?2,  decl.  = 
— 19°26'28"  (J2000.0)  (just  west  of  the  SNR  shell  structure) 
with  a  positive  spectral  index  has  also  been  detected  in  a  3  cm 
radio  recombination  line  (RRL)  survey  by  Lockman  (1989) 
with  the  NRAO  140  foot  (43  m)  telescope  with  a  3'  beam. 
This  RRL  has  a  velocity  of  18.5  km  s-1  and  Av  ~  24  km  s_1. 
Other  tracers  of  star  formation  have  also  been  detected  in  this 
region  with  similar  velocities,  including  a  molecular  cloud  in 
CS  and  methanol  masers  (Bronfman,  Nyman,  &  May  1996; 
Walsh  et  al.  1998).  Confusion  from  the  associated  H  n  region 
likely  explains  why  this  SNR  has  not  been  discovered  in 
lower  resolution  surveys. 

3.1.5.  New  SNR  C:  Gil.  18+0.11 

Figure  6  shows  the  integrated  continuum  spectrum  for  SNR 
C  at  G1 1.18+0. 11,  while  Figure  le  shows  the  morphology 
of  the  spatially  resolved  330/1465  MHz  spectral  index. 
The  integrated  spectral  index  for  this  SNR  candidate  is 
—0.44  ±  0.01.  The  region  of  positive  spectral  index  near  the 
northern  end  of  candidate  C  (J2000.0;  R.A.  =  18h09m54s, 
decl.  =  —19°  10'  12")  is  coincident  with  a  RRL  detected  by 
Lockman  (1989),  with  a  velocity  of  7.3  km  s-1  and  Ad  ^54  km 
s-1  using  the  NRAO  140  foot  at  3  cm  (to  within  their  3'  beam). 
It  is  difficult  to  estimate  the  contribution  of  this  thermal 
component  to  the  integrated  SNR  flux  because  of  the  absence 
of  spatially  resolved  data  at  a  higher  frequency.  For  example, 
an  unresolved  source  is  apparent  at  the  position  of  the  SNR 
candidate  in  the  2695  MHz  Bonn  survey  images,  but  it  is 
confused  with  both  the  H  n  region  traced  by  the  RRL,  as  well  as 
the  steep-spectrum  point  source  just  east  of  the  SNR  shell  that 
is  likely  extragalactic  (see  Fig.  le).  For  this  reason,  no  effort 
was  made  to  remove  this  thermal  component  from  the 
integrated  flux  densities,  and  the  reported  a  should  be  viewed 
as  a  lower  limit. 

4.  DISCUSSION 

4.1.  The  Gil. 2— 0.3  Field  at  Other  Wavelengths 

Comparisons  of  high-resolution  radio  and  infrared  Galactic 
plane  images  are  quite  useful  in  distinguishing  different 
components  of  the  ISM  (see,  for  example,  Cohen  &  Green 
2001).  Figure  8  shows  a  ~20"  resolution  mid-infrared  8.28  /mi 
image  from  the  Midcourse  Space  Experiment  Satellite  ( MSX ) 
Galactic  Plane  Survey  of  the  Gil. 2— 0.3  field  with  the  4  mJy 


beam-1  1465  MHz  contour  superposed  (see  Price  et  al.  2001 
for  details  about  the  MSX  image  processing).  In  the  Galactic 
plane  the  8.28  pm  emission  is  dominated  by  contributions 
from  polycyclic  aromatic  hydrocarbon  molecular  line  emis¬ 
sion  and  by  thermal  emission  from  warm  dust  around  evolved 
stars  and  regions  of  star  formation  (normal,  and  ultracompact 
H  ii  regions).  The  details  of  the  origin  of  the  diffuse  mid-IR 
emission  is  less  certain  but  may  be  associated  with  the  dif¬ 
fuse  warm  ionized  medium  (WIM).  The  1465  MHz  radio 
continuum  image  has  contributions  from  both  nonthermal, 
as  well  as  ionized  thermal  sources.  In  Figure  8  sources  with 
both  strong  infrared  and  radio  detections  are  likely  to  be 
normal  or  ultracompact  H  n  regions;  unresolved  sources  with 
radio  but  no  infrared  emission  are  nonthermal  and  are  most 
likely  extragalactic  but  could  also  be  very  young  SNRs, 
pulsars,  or  unresolved  PWNe;  unresolved  sources  with 
infrared,  but  no  radio  continuum,  emission  are  most  likely 
normal  and  evolved  stars;  extended  infrared  sources  without 
radio  continuum  emission  trace  the  very  youngest  stages  of 
star  formation  where  the  surrounding  medium  has  not  yet  been 
ionized  (regions  that  are  optically  thick  even  at  8.28  pm  can 
appear  as  mid-IR  voids  with  no  radio  emission,  e.g.,  in  the 
southwest  comer  of  Fig.  8);  and,  finally,  extended  radio  con¬ 
tinuum  emission  without  strong  infrared  emission  must  be 
nonthermal  and  may  be  SNRs  (see  Cohen  &  Green  2001,  for  a 
more  detailed  discussion  of  these  classifications  and  some 
caveats).  There  are  a  few  cases  where  extended  infrared  emis¬ 
sion  has  been  detected  toward  SNRs,  but  this  only  seems  to 
occur  when  an  SNR  shocks  a  nearby  molecular  cloud  (i.e.,  3C 
391,  Reach  et  al.  2002). 

A  low  ratio  of  infrared-to-radio  continuum  flux  has  been 
used  for  some  time  as  a  tool  to  identify  SNR  candidates  in 
Galactic  plane  surveys  (Whiteoak  &  Green  1996).  However, 
at  the  comparatively  poorer  resolutions  used  in  past  surveys 
(>1')  application  of  this  test  could  be  ambiguous  in  confused 
regions  of  both  thermal  and  nonthermal  emission.  The  high- 
resolution  images  shown  in  Figure  8  dramatically  demonstrate 
the  location  of  the  Galactic  plane  running  through  the  eastern 
half  of  the  image,  as  well  as  the  confusion  associated  with 
star  formation  in  this  region,  particularly  near  the  new  SNRs 
B:  Gil. 03— 0.05  and  C:  Gll.18+0.11.  Many  of  the  sources  in 
our  330/1465  MHz  spectral  index  maps  (Figs,  la-le)  with 
a  >  0.1  have  counterparts  at  8.28  pm  and  are  likely  to  be  H  n 
regions.  The  radio  point  sources  located  toward  the  southwest 
comer  of  G1 1.4— 0.1  and  to  the  east  of  SNR  C  (see  §§  3.1.2 
and  3.1.5)  are  not  coincident  with  infrared  emission  and, 
therefore,  are  statistically  likely  to  be  extragalactic  although  as 
mentioned  above  they  could  also  be  pulsars.  No  extended 
infrared  features  (above  the  diffuse  Galactic  infrared  back¬ 
ground)  are  detected  toward  the  five  SNRs  discussed  here. 
There  is  also  no  evidence  for  an  infrared  counterpart  to  the 
extended  linear  radio  continuum  feature  that  lies  tangent  to  the 
southern  edge  of  Gil. 2— 0.3  (§  3.1.1). 

The  ASCA  and  ROSATA-xay  archives  were  also  searched  for 
possible  correlations  with  the  SNRs  in  the  G1 1.2— 0.3  field. 
With  the  exception  of  Gil. 2— 0.3  itself,  whose  X-ray  proper¬ 
ties  have  been  discussed  by  a  number  of  authors  (see,  e.g., 
Vasisht  et  al.  1996),  none  of  the  other  SNRs  in  the  field  were 
detected  by  either  ASCA  or  ROSAT.  However,  there  is  a  region 
of  diffuse  X-ray  emission  in  between  SNRs  B:  G1 1.03—0.05 
and  C:  G1 1.18+0.1 1  in  an  image  from  the  ASCA  Galactic  plane 
survey  (47  ks  exposure  with  the  GIS  instrument;  also  see 
Bamba  et  al.  [2003]).  An  ASCA  X-ray  image  spanning  the  en¬ 
ergy  range  2-10  keV  in  the  region  of  SNRs  B  and  C  is  shown 
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Fig.  8. — MSX  8.28  //m  image  of  the  G1 1.2— 0.3  field  with  the  4  mjy  beam-1  1465  MHz  contour  superposed.  The  resolution  of  the  8.28  //m  image  is  ~20".  Note 
that  many  of  the  sources  with  positive  spectral  indices  in  Fig.  7  are  coincident  with  thennal  gas  traced  by  strong  8.28  fi m  emission. 


in  Figure  9.  The  X-ray  data  were  processed  using  the  technique 
described  in  Roberts,  Romani,  &  Kawai  (2001).  The  origin  of 
the  diffuse  X-ray  emission  near  the  center  of  Figure  9  is  un¬ 
certain,  however,  the  radio  pulsar  PSR  J1809— 1917  described 
in  §  4.3  is  coincident  with  a  peak  in  the  X-ray  emission  (see 
Fig.  9).  Low  signal-to-noise  diffuse  emission  is  also  coincident 
with  the  middle  and  northern  parts  of  G1 1 .03—0.05.  The  X-ray 
spectra  derived  from  the  ASCA  data  lack  sufficient  counts 
to  adequately  fit  their  spectra.  It  would  be  interesting  to 
investigate  this  region  in  the  X-ray  with  higher  sensitivity  to 
determine  the  nature  of  the  diffuse  emission. 

4.2.  Distances  from  the  X-D  relation 

The  Tj-D  relation  (surface  brightness  vs.  diameter)  has  long 
been  used  to  assign  distances  to  SNRs  without  alternative 
distance  estimates,  based  on  the  observed  properties  of  SNRs 
at  known  distances.  Such  relations  rely  on  the  assumption  that 
supernova  explosion  energies,  ambient  densities,  and  the 
strength  and  evolutionary  state  of  their  magnetic  fields  are 
similar.  It  remains  unclear  to  what  extent  these  assumptions  are 
true,  particularly  in  regard  to  the  ambient  density.  Indeed, 
while  X-D  relations  often  provide  a  fairly  tight  correlation  for  a 
statistical  sample,  their  application  to  individual  SNRs  is  fairly 
dubious  (see  Case  &  Bhattacharya  1998,  for  further  dis¬ 
cussion).  Despite  this  caveat,  we  have  used  the  most  recent 


determination  of  the  Galactic  X-D  relation  derived  by  Case  & 
Bhattacharya  (1998): 

X1GHz  =  2.07±f'i2  x  10 -17Z)(-2-38±0-26),  (2) 

to  estimate  the  distances  to  the  five  SNRs  discussed  in  this 
paper  (see  Table  2).  The  sample  used  to  derive  the  fit  in 
equation  (2)  includes  data  from  36  Galactic  shell-type  SNRs. 
Case  &  Bhattacharya  (1998)  find  that  the  average  fractional 
deviation  of  an  individual  SNRs  distance  in  their  sample  from 
that  calculated  from  equation  (2)  is  ~40%. 

Table  2  shows  the  parameters  used  to  derive  the  SNR 
distances.  The  distance  to  G1 1.2— 0.3  estimated  from  equation 
(2),  6  kpc,  is  fairly  close  to  that  estimated  from  H  i  absorption 
(see  §  3.1.1).  Note  that  Gil. 2— 0.3  was  not  included  in  the 
Case  &  Bhattacharya  derivation  of  the  X-Z)  relation.  The 
estimated  distances  to  the  other  SNRs  are  G1 1.4— 0.1:  9  kpc, 
Gil. 15— 0.71:  24  kpc,  Gll.03-0.05:  16  kpc,  and  Gil. 18+ 
0.11:  17  kpc.  The  range  of  these  values,  assuming  a  40% 
fractional  error,  is  also  listed  in  Table  2.  It  is  notable  that  the 
fraction  of  low  surface  brightness  compared  with  high  surface 
brightness  shell-type  remnants  among  the  known  SNRs  is 
much  higher  than  the  fraction  used  to  derive  equation  (2)  (i.e., 
few  SNRs  with  low  surface  brightness  have  known  distances; 
Case  &  Bhattacharya  1998).  Hence,  equation  (2)  is  significantly 
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Fig.  9. — ASCA  X-ray  image  in  the  2-10  keV  range  of  the  region 
containing  the  SNRs  B:  G1 1.03—0.05  and  C:  G1 1.18+0.11  (see  Figs.  1  a  and 
1  b)  with  the  4,  10,  and  20  mJy  beam-1  1465  MHz  contours  superposed.  The 
X-ray  data  have  been  convolved  to  the  approximate  FWHM  resolution  of  the 
GIS  detector  of  50".  The  white  asterisk  shows  the  location  of  pulsar  PSR 
J 1809— 1917.  The  image  was  processed  as  described  in  Roberts  et  al.  2001. 

biased  toward  high  surface  brightness  shell-type  SNRs.  The 
effect  of  this  bias  on  equation  (2),  and  hence  our  three  new  low 
surface  brightness  shell-type  remnants  are  uncertain.  It  is 
interesting  that  the  high  end  of  the  distance  ranges  presented  in 
Table  2  (assuming  a  40%  fractional  error),  would  unrealistically 
place  all  three  of  the  new  SNRs  outside  the  bounds  of  the 
Galaxy,  perhaps  suggesting  that  the  distance  to  these  SNRs  are 
not  well  fitted  by  equation  (2).  Also,  note  that  the  cases  of 
Gil. 15— 0.71  and  Gil. 18+0. 11  are  particularly  uncertain 
(independent  of  the  accuracy  of  eq.  [2])  since  the  angular  size 
( 6 )  of  these  two  SNRs  are  difficult  to  estimate  since  they  only 
appear  as  partial  shells  in  the  current  images.  This  analysis 
suggests  that  the  three  new  SNRs  may  lie  on  the  far  side  of  the 
Galaxy,  but  H  i  absorption  measurements  are  needed  to  better 
constrain  their  distances. 

4.3.  Pulsar  Association? 

In  addition  to  the  pulsar  associated  with  Gil. 2— 0.3  de¬ 
scribed  in  §  3.1.1,  one  other  pulsar  has  been  discovered  within 
the  confines  our  low-frequency  field  of  view.  PSR  J1809— 1917 
was  discovered  in  the  Parkes  Multibeam  Pulsar  Survey  by 


Morris  et  al.  (2002)  and  is  located  at  R.A.  =  18h09m43®2, 
deck  =  — 19°17'38"  (J2000.0).  The  position  of  the  pulsar  is 
indicated  on  Figure  9.  It  has  a  pulse  period  of  83  ms,  a 
flux  density  at  1400  MHz  of  2.5  mJy,  dispersion  measure 
DM  =  197  pc  cm-3,  a  dispersion  measure  distance  of  —  4  kpc, 
and  a  characteristic  age  of  50  kyr.  We  do  not  detect  the  pulsar 
at  either  330  or  1465  MHz.  This  is  not  surprising  for  the 
1465  MHz  data  because  the  noise  rises  to  —2  mJy  beam-1  near 
the  edge  of  the  image  (due  to  the  primary  beam  correction). 
However,  if  we  assume  a  typical  pulsar  spectral  index  of  —2  the 
predicted  330  MHz  flux  density  is  45  mJy  beam-1  which  should 
be  easily  detectable  in  our  330  MHz  image  (rms  noise  3  mJy 
beam-1).  Alternatively,  the  pulsar’s  spectral  index  could  be 
steeper  than  average  (a  =  —  1  predicts  S330  =  11  mJy  beam-1). 
In  the  current  ATNF  pulsar  catalog,11  13%  of  the  283  pulsars 
with  a  measured  spectral  index  have  a  >  —  1 . 

The  characteristic  age  estimate  (50  kyr)  assumes  that  the 
birth  spin  period  is  negligible  compared  with  the  current  spin 
period.  However,  recent  studies  on  a  small  number  of  pulsars 
for  which  an  independent  age  estimate  can  be  made  suggest 
that  the  birth  period  is  actually  between  20  and  130  ms  (see, 
e.g.,  Kaspi  et  al.  2001;  Kramer  et  al.  2003).  Thus,  for  young 
pulsars  that  currently  have  short  spin  periods  (i.e.,  little 
spindown  has  yet  occurred),  a  fundamental  assumption  of  the 
characteristic  age  calculation  is  violated.  This  implies  that  the 
true  ages  of  these  pulsars  could  be  significantly  younger  than 
suggested  by  the  characteristic  age.  For  example,  if  the  birth 
period  of  PSR  J 1809— 1917  was  —50  ms  and  we  assume  that 
the  spindown  is  dominated  by  dipole  radiation  with  a  braking 
index  of  3,  then  the  true  age  of  the  pulsar  is  ~25  kyr,  half  the 
characteristic  age  (see  Kramer  et  al.  2003). 

PSR  J 1 809 — 1917  is  located  about  8.' 5  from  the  center 
of  Gil. 03— 0.05  (new  SNR  B)  and  6'  from  the  center  of 
G1 1.1 8+0.1 1  (new  SNR  C).  If  these  SNRs  are  located  —  17  kpc 
away  (see  Table  2),  compared  with  the  —  4  kpc  distance 
estimated  for  the  pulsar,  then  obviously  there  can  be  no 
association.  However,  as  described  in  §  4.2,  the  uncertainty  in 
distances  derived  from  the  H-D  relation  is  considerable.  Also, 
given  that  the  pulsar  is  fairly  young  and  the  notorious 
difficulty  with  establishing  SNR/pulsar  associations,  it  is 
worthwhile  to  determine  whether  the  pulsar  could  have 
reasonably  moved  to  its  current  location  from  the  center  of 
either  remnant  (assuming  they  are  at  the  same  distance). 
No  independent  distance  has  yet  been  obtained  for  PSR 
J1 809— 1917,  so  we  will  assume  the  dispersion  measure  pul¬ 
sar  distance  of  4  kpc.  This  distance  is  based  on  a  model  for 
the  electron  distribution  along  the  line  of  sight  to  the  pulsar 
and  is  probably  good  to  within  —25%  (Taylor  &  Cordes 
1993).  Assuming  that  the  pulsar  is  50  kyr  old,  the  pulsar’s 
transverse  speed  would  need  to  be  —200  km  s-1  to  originate 
from  Gil. 03— 0.05  and  -140  km  s-1  for  Gil. 18+0.11.  Such 
transverse  velocities  are  quite  reasonable  for  pulsars  and 
woidd  still  be  plausible  if  the  true  age  of  the  pulsar  is  younger 
than  the  characteristic  age  as  described  above  or  the  pulsar  is 
more  distant.  For  example,  from  analysis  of  the  velocity 
distributions  of  a  large  sample  of  pidsars  Arzoumanian, 
Chemoff,  &  Cordes  (2002)  predict  that  —50%  of  pulsars  have 
velocities  ^500  km  s-1.  Further  assessment  of  the  reality  of  a 
possible  association  with  either  remnant  will  have  to  await 
more  reliable  distance  estimates  for  the  SNRs. 


11  See  http://www.atnf.csiro.au/researchypulsar/psrcat. 
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TABLE  2 

Derived  SNR  Parameters 


■Si ghz  9  si  GHz  Distance15  Range0 

Source  Name  aa  (Jy)  (arcmin)  (W  m“2  Hz  Sr  ')  (kpc )  (kpc) 


G1 1.2— 0.3 .  -0.57  19.3d  4.5  1.4  x  10“19  6  4-10 

Gll.4-0.1 .  -0.59  7.0  7  x  9  1.7  x  Id  •’:l  9  6-14 

Gil. 15-0.71  (A)....  -0.82  0.6  8  x  7  1.7  X  KT21  24  17-40 

Gll.03-0.05  (B)....  -0.60  1.6  6  x  7  5.9  X  10“21  16  12-27 

Gil. 18+0.11  (C) .  -0.44  1.5  8  x  5  5.8  X  10“21  17  12-28 


a  Fitted  spectral  indices  from  integrated  flux  measurements  shown  in  Figs.  2-6. 
b  Distances  calculated  using  eq.  (2). 

0  Range  of  distances  assuming  a  40%  fractional  error. 
d  After  subtracting  1  Jy  for  the  PWN. 


4.4.  Low-Frequency  Turnovers  for  Gll.2—0.3  and  Gll.4—0.1 

The  integrated  continuum  spectra  shown  in  Figures  2  and 
3  for  Gll.2—0.3  and  Gll.4—0.1  show  clear  signs  of  a  low- 
frequency  turnover.  Their  optical  depths  at  74  MHz  of  0.56 
and  1.1,  respectively,  are  higher  than  those  found  by  Kassirn 
(1989a)  by  about  40%.  Because  our  new  spectra  include 
critical  new  data  points  at  74  and  330  MHz  and  we  have 
chosen  to  only  use  integrated  flux  densities  from  the  lit¬ 
erature  with  error  estimates  that  are  less  than  20%,  our  esti¬ 
mated  optical  depths  are  more  accurate.  The  three  new  SNRs 
also  show  indications  of  a  low-frequency  turnover  based  on 
their  74  MHz  upper  limits,  which  we  hope  to  confirm  with 
future  higher  resolution  74  MHz  data.  From  a  survey  of  the 
integrated  continuum  spectra  of  47  Galactic  SNRs,  Kassim 
(1989a)  found  that  about  two-thirds  of  the  observed  rem¬ 
nants  show  a  low-frequency  turnover  below  Si  200  MHz. 
This  turnover  is  thought  to  be  due  to  free-free  absorption  by 
intermediate  temperature  (~3000-8000  K)  and  density  (~1- 
10  cnr3)  ionized  thermal  gas  along  the  line  of  sight  to  the 
SNRs  (see,  e.g.,  Kassim  1989a  and  references  therein). 
Based  on  the  fact  that  such  absorption  is  not  seen  toward 
all  SNRs  and  does  not  seem  to  depend  on  distance,  it  is 
unlikely  that  the  absorbing  medium  is  a  distributed  com¬ 
ponent  of  the  ISM  and  is  more  likely  associated  with  discrete 
sources.  Kassim  estimates  that  the  filling  factor  of  ionized 
material  with  the  required  electron  densities  and  temperatures 
is  <1%. 

A  number  of  origins  for  the  absorbing  medium  have  been 
suggested.  Two  obvious  sources  of  Galactic  ionized  gas,  the 
warm  ionized  medium  (WIM)  and  H  n  regions,  do  not  provide 
the  appropriate  physical  conditions.  The  WIM  has  electron 
densities  that  are  too  low  by  a  factor  of  ~10  to  100,  and  its 
filling  factor  is  too  high  (>10%;  Kulkami  &  Heiles  1987). 
Normal  H  n  regions  have  electron  densities  and  temperatures 
that  are  far  too  high  to  account  for  the  moderate  low-fre¬ 
quency  optical  depths  that  are  observed,  and  they  would  be 
easily  detected  in  higher  frequency  continuum  observations. 
One  reasonable  suggestion  is  that  the  extended  envelopes  of 
normal  H  n  regions  (EHEs)  might  well  provide  the  requisite 
temperatures  and  densities  (Anantharamaiah  1985,  1986).  It  is 
unclear  whether  such  envelopes  extend  far  enough  into  the 
ISM  to  account  for  a  large  fraction  of  the  SNRs  showing  low- 
frequency  absorption. 

Alternatively,  Heiles,  Reach,  &  Koo  (1996)  suggest  that 
Galactic  “worms”  may  cause  absorption  toward  some  sources. 
Worms  are  thought  to  be  created  after  the  first  generations  of 


an  OB  star  cluster  explode  as  supemovae,  blowing  a  large 
evacuated  superbubble,  the  walls  of  which  are  then  ionized 
by  later  generations  of  OB  stars.  The  spatial  extent  of  such 
worms  can  be  several  degrees  across.  Koo,  Heiles,  &  Reach 
( 1 992)  have  cataloged  1 1 8  Galactic  worm  candidates  from  single¬ 
dish  H  i  and  radio  continuum  data  plus  correlations  with 
IRAS  infrared  data.  All  of  the  worm  candidates  lie  toward  the 
inner  Galaxy.  Heiles  et  al.  (1996)  find  that,  in  some  cases, 
worms  also  emit  RRLs  and  that  the  electron  densities  and 
temperatures  expected  for  worms  are  commensurate  with 
those  needed  to  produce  free-free  absorption  at  low  radio 
frequencies. 

The  low-frequency  absorption  observed  toward  the  SNR 
W49  B  (G43.3— 0.2)  by  Lacey  et  al.  (2001)  could  be  a  case  of 
absorption  by  a  worm.  An  early  identification  of  the  existence 
of  ionized  thermal  gas  toward  W49  B  came  from  the  detection 
of  RRL  lines  at  ^60  km  s_1  (Pankonin  &  Downes  1976). 
Lacey  et  al.  (2001)  subsequently  found  excellent  agreement 
between  the  spatially  resolved  morphology  of  the  free-free 
absorption  at  74  MHz  and  that  ofH  i  absorption  at  ~60  km  s  1 . 
There  are  no  obvious  H  ii  regions  near  W49B  (in  the  plane  of 
the  sky)  with  a  velocity  of  ~60  km  s_1  that  could  be  candidates 
for  producing  an  EHE.  Lacey  et  al.  suggest  that  the  absorption  is 
due  to  a  few  compact  H  n  regions  (and  their  envelopes)  along 
the  line  of  sight  to  W49  B  itself,  but  it  is  difficult  to  see  how  the 
cores  of  such  H  n  regions  could  escape  detection  at  higher 
frequencies  (Moffett  &  Reynolds  1994).  Interestingly,  one  of 
the  worms  cataloged  by  Koo  et  al.  (1992)  is  in  this  general 
direction  (GW  44.8—1.8),  has  an  extent  of  ~2°,  and  has  ve¬ 
locities  between  43  to  66  km  s~ 1 .  Future  RRL  observations  may 
determine  whether  the  physical  conditions  in  the  worm  are  in 
agreement  with  those  needed  to  produce  the  free-free 
absorption. 

It  is  also  possible  that  some  SNRs  could  create  their  own 
absorbing  medium  if  they  happen  to  run  into  a  nearby  molec¬ 
ular  cloud,  assuming  that  the  shock  is  a  fast,  ionizing  J-type 
shock.  Evidence  for  this  case  has  recently  been  discovered  for 
the  SNR  3C  391,  which  is  known  to  be  interacting  with  a 
nearby  molecular  cloud  (see,  e.g.,  Reach  &  Rho  1999).  The 
spatially  resolved  morphology  of  the  74  MHz  free-free 
absorption  toward  3C  391  matches  very  well  with  that  of 
mid-infrared  fine  structure  atomic  lines,  which  are  known  to 
be  good  tracers  of  ionic  shocks  (Brogan  et  al.  2004;  Reach 
et  al.  2002).  Since  all  of  these  phenomena  exist,  it  is  likely  that 
each  of  them  is  responsible  for  at  least  some  fraction  of  the 
observed  free-free  absorption,  and  only  by  studying  a  number 
of  sources  in  detail  can  the  dominant  process  be  determined. 
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The  emission  measures  implied  by  our  newly  derived  op¬ 
tical  depths  for  Gil. 2— 0.3  and  Gil. 4— 0.1  can  be  calcu¬ 
lated  from 

t{v)  =  1.643  x  \tfa(Te,  v)v~1A  EM  Tf135,  (3) 

where  EM  is  the  emission  measure,  Te  is  the  electron 
temperature,  and  a(Te,  v )  is  the  Gaunt  factor  and  is  ~  1  at  the 
temperatures  and  densities  discussed  here  (Dulk  &  Slee  1975). 
Assuming  an  electron  temperature  of  5000  K,  the  average 
emission  measure  (EM)  for  absorbing  gas  toward  G1 1.2— 0.3 
and  Gil. 4— 0.1  at  74  MHz  is  2.8  xlO3  cm  6  pc  and 
5.6  xlO3  cm~6  pc,  respectively.  Accounting  for  the  lower 
Te  used  here,  5000  K  instead  of  10,000  K,  these  emission 
measures  are  similar  to  those  found  toward  the  SNR  W49  B 
(Lacey  et  al.  2001). 

There  are  no  indications  in  the  spectral  index  maps  of 
Gil. 2-0.3  and  G1 1.4-0. 1  between  330  and  1465  MHz  of  the 
morphology  of  the  absorbing  medium.  This  is  not  surprising 
given  that  the  330  MHz  optical  depths  are  only  0.02  and  0.05 
for  Gil. 2— 0.3  and  G1 1.4— 0.1,  respectively.  There  is  some 
indication  in  our  145"  x  92"  resolution  74  MHz  image  that 
Gil. 2— 0.3  suffers  from  more  absorption  on  its  eastern  side 
than  the  western  side,  but  this  will  need  to  be  confirmed  by 
higher  resolution  74  MHz  data.  The  H  n  region  Gil.  11— 0.4 
located  just  south  of  Gil. 2— 0.3  (see  Fig.  la),  is  a  potential 
candidate  to  provide  an  absorbing  EHE  if  its  distance  along 
the  line  of  sight  is  closer  to  us  than  that  of  the  remnant.  The 
kinematics  of  the  molecular  gas  associated  with  this  H  n 
region  are  consistent  with  a  location  in  the  3  kpc  arm  of  the 
Galaxy,  suggesting  that  the  distance  to  the  H  n  region  is 
~5  kpc,  similar  to  that  of  the  SNR  (Solomon  et  al.  1987; 
Kurtz,  Churchwell,  &  Wood  1994).  A  detailed  comparison  of 
H  i  absorption  toward  both  Gil. 2— 0.3  and  Gil.  11— 0.4  will 
be  needed  to  better  determine  their  relative  distances.  No 
known  H  n  regions  are  in  close  enough  proximity  (in  the  plane 
of  the  sky)  to  Gil. 4— 0.1  (see  Fig.  la)  to  provide  an  EHE 
candidate  for  this  SNR. 

Spatially  resolved  74  MHz  images  have  been  crucial  in 
determining  the  nature  of  the  absorbing  medium  for  W49  B 
(Lacey  et  al.  2001)  and  3C  391  (Brogan  et  al.  2004).  Future 
spatially  resolved  74  MHz  images  are  likely  to  be  equally 
important  in  identifying  the  source  of  ionized  thermal  gas 
toward  Gll.2-0.3  and  Gll.4-0.1. 

4.5.  Implications  for  the  Number  of  Galactic  SNRs 

A  number  of  recent  studies,  ranging  from  analyzing  the  rate 
of  historical  Galactic  supemovae,  the  rate  of  supemovae  in  the 
Local  Group  of  galaxies,  to  the  birth  rate  of  OB  stars  within 
1  kpc  of  the  Sun,  suggest  that  the  rate  of  Galactic  supemovae 
is  about  one  per  50  yr  (see,  e.g.,  Tammann  et  al.  1994).  If  we 
assume  that  supernova  remnants  can  be  distinguished  from  the 
Galactic  background  nonthermal  emission  for  about  50,000  to 
100,000  yr,  then  there  should  be  on  the  order  of  1000  to  2000 
SNRs  in  our  Galaxy  today.  This  number  is  far  above  the 
number  of  currently  known  Galactic  SNRs,  ~230  (Green 
2002). 

For  some  time  observational  selection  effects  have  been 
thought  to  be  responsible  for  the  low  number  of  Galactic 
SNRs  currently  known  (see,  e.g.,  Green  1991).  These  selection 
effects  include  (1)  the  difficulties  associated  with  detecting  the 
more  mature,  large,  faint  remnants  in  the  inner  Galaxy  due  to 
confusion  with  the  Galactic  synchrotron  background  radiation; 
(2)  the  challenges  involved  with  identifying  small,  young 


shell-type  SNRs  that  can  be  easily  confused  with  steep  spec¬ 
trum  extragalactic  sources;  (3)  difficulties  with  distinguishing 
young  and  older  shell-type  remnants  when  they  are  confused 
with  H  n  regions  that  lie  along  the  line  of  sight;  and  (4)  the  fact 
that  about  5%  of  SNRs,  known  as  plerions,  have  no  identifiable 
shell  and  have  flat  spectra,  so  that  they  are  difficult  to  dis¬ 
tinguish  from  H  ii  regions  even  in  high-resolution  images  (see, 
e.g.,  Helfand  et  al.  1989).  The  solution  to  the  first  problem  is 
obviously  to  make  high  surface  brightness  sensitivity  obser¬ 
vations,  while  the  second  requires  high  spatial  resolution  in  the 
hopes  of  resolving  a  candidate  SNRs  shell  morphology.  The 
solution  to  number  (3)  is  to  obtain  equally  sensitive,  high- 
resolution  multiwavelength  data  so  that  sources  of  thermal 
emission  can  be  separated  from  nonthermal  emission.  The  so¬ 
lutions  to  number  (4)  include  looking  for  strong  linear  polar¬ 
ization  or  an  associated  pulsar. 

In  the  ~1  deg2  field  centered  on  G1 1.2— 0.3  we  have 
overcome  the  first  three  of  the  above-mentioned  problems  by 
obtaining  high  surface  brightness  sensitivity,  high-resolution, 
multiwavelength  data.  These  data  have  resulted  in  the  dis¬ 
covery  of  three  new  SNRs.  All  three  new  remnants  appear 
fairly  old  (i.e.,  they  are  relatively  extended  and  have  low  sur¬ 
face  brightnesses)  and,  based  on  the  somewhat  dubious  Y,-D 
relation,  are  located  on  the  other  side  of  the  Galaxy  (16  to 
24  kpc).  We  have  not  identified  any  new  young  SNRs  in  this 
field.  However,  statistically  this  is  not  surprising  since  this 
field  also  contains  the  historical  remnant  Gil. 2— 0.3  from 
A.D.  386,  and  we  have  only  sampled  a  small  fraction  of  the 
inner  Galactic  plane.  If  we  define  “young”  SNRs  as  those 
having  ages  less  than  2000  yr,  then  assuming  one  SNR  every 
50  yr  implies  that  there  are  ^40  such  “young”  remnants  in  the 
Galaxy.  There  are  only  six  known  historical  remnants 
(meaning  that  sightings  of  their  optical  emission  were 
recorded  in  historical  documents):  the  Crab,  SN1006,  3C  58, 
Tycho,  Kepler,  and  Gil. 2— 0.3,  plus  a  handful  of  suspected 
“young”  SNRs  including  CasA.  Thus,  there  are  on  the  order 
of  30  “young”  SNRs  that  have  yet  to  be  discovered.  These 
SNRs  are  presumably  missing  from  the  historical  records 
because  the  optical  emission  from  their  initial  explosions  was 
obscured  by  intervening  material.  Indeed,  five  of  the  six  his¬ 
torical  SNRs  described  above  are  a  degree  or  more  off  the 
plane. 

Extrapolating  from  this  approximately  1  deg2  field  of 
view,  we  can  estimate  how  many  SNRs  should  be  detected  in 
our  planned  330  MHz  survey  of  the  inner  Galaxy,  from 
+4°  <  /  <  +20°  with  \b\  <  1°.  For  this  estimate,  we  make 
the  simplifying  assumption  that  the  distribution  of  SNRs  within 
this  strip  of  the  Galactic  plane  is  fairly  uniform.  However, 
this  assumption  could  be  violated  if  some  lines  of  sight  pass 
through  more  spiral  arms  than  others.  Cordes  &  Lazio  (2003) 
have  recently  completed  a  new  analysis  of  the  Galactic  free 
electron  distribution  determined  from  pulsar  dispersion 
measurements  (this  work  extends  and  improves  upon  the 
original  study  by  Taylor  &  Cordes  1993).  From  the  Galactic 
electron  distribution  and  many  other  lines  of  evidence,  these 
authors  have  also  developed  an  updated  model  for  the  Galac¬ 
tic  spiral  arm  structure.  The  Cordes  &  Lazio  (2003)  spiral 
arm  model  suggests  that  there  are  two  potential  disparities  in 
the  total  number  of  spiral  arms  in  the  longitude  range 
+4°  <  /  <  +20°.  These  include  the  fact  that  (1)  the  Norma 
arm  does  not  appear  to  extend  beyond  /  >  10°  on  the  near  side 
of  the  Galactic  center;  and  (2)  the  Crux-Scutum  arm  does  not 
appear  to  extend  to  longitudes  smaller  than  l  =  18°  on  the  far 
side  of  the  Galactic  center.  An  even  greater  discrepancy  is 
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possible  if  a  line  of  sight  passes  through  a  tangent  point  of  a 
spiral  arm  where  gas,  H  n  regions,  and,  presumably,  SNRs 
tend  to  pile  up.  This  may  well  be  true  for  lines  of  sight  near 
/  =  20°.  In  any  case;  the  region  near  I  =  11°  2  does  not  appear 
to  be  a  particularly  special  one  within  the  +4°  <  l  <  +20° 
range,  suggesting  that  this  field  of  view  should  not  be 
overdense  with  SNRs  compared  with  the  rest  of  the  survey 
region. 

The  above  caveats  aside,  the  results  from  our  ~  1  deg2  field 
toward  Gil. 2— 0.3  suggest  that,  near  b  —  0°,  the  density  of 
SNRs  in  this  part  of  the  Galactic  plane  could  be  as  high  as  five 
per  square  degree.  We  will  sample  approximately  32  deg2  in 
our  survey,  but  we  assume  that  only  the  ~16  deg2  near  b  =  0° 
(i.e.,  \b\  ;S  0?5)  has  an  SNR  density  of  five  per  square  degree, 
with  the  extra  area  accounting  somewhat  for  density  varia¬ 
tions.  These  assumptions  imply  a  total  SNR  count  in  the 
+4°  %  l  %  +20°  and  \b\  5=  1°  region  of  ~80.  This  number 
is  in  reasonable  agreement  with  that  predicted  by  Helfand  et  al. 
(1989)  for  this  region  of  the  Galactic  plane  (63)  based  on  a 
Galaxy-wide  treatment  of  SNR  densities  (see  also  Case  & 
Bhattacharya  1998  and  references  therein).  Within  the  bounds 
of  our  future  survey  region  only  14  SNRs  are  currently  known 
(Green  2002)!  Only  the  actual  results  from  our  more  extensive 
survey  will  determine  whether  this  estimate  is  realistic,  but,  in 
any  case,  it  will  go  a  long  way  toward  determining  the  true 
distribution  of  Galactic  SNRs. 

5.  CONCLUSIONS 

In  a  ~  1  deg2  field  of  view  centered  on  the  known  Galactic 
SNR  G1 1.2— 0.3,  we  have  identified  three  new  SNRs.  Pre¬ 
viously,  this  field  was  thought  to  contain  only  two  SNRs.  The 
integrated  spectral  indices  of  the  new  SNRs:  A:  G1 1.15— 0.71, 
B:  Gil. 03— 0.05,  and  C:  Gil. 18+0.11  range  from  -0.80 
to  —0.44.  Two  of  the  new  SNRs  (G1 1.03—0.05  and 
Gil.  18+0.11)  have  confusing  H  n  regions  nearby  (along  the 
line  of  sight)  that  are  not  spatially  resolved  from  the  SNRs  at 
the  resolution  of  the  Bonn  2695  MHz  survey,  and  this  fact  is 
likely  responsible  for  their  not  having  been  previously 
identified.  We  have  also  confirmed  the  low-frequency  spectral 
turnover  in  the  integrated  radio  continuum  spectra  of 
Gil. 2— 0.3  and  G1 1.4— 0.1  previously  observed  by  Kassim 
(1989a)  that  is  due  to  free-free  absorption  along  the  line  of 
sight.  The  addition  of  our  new  data  suggest  that  the  optical 
depths  toward  these  SNRs  at  74  MHz  are  0.56  for  Gil. 2— 0.3 
and  1.1  for  G1 1.4— 0.1.  The  upper  limits  for  the  74  MHz 
integrated  flux  densities  for  the  three  new  SNRs  also  suggest 


that  these  remnants  suffer  from  absorption.  While  quite 
uncertain,  a  recent  formulation  of  the  £-£>  relation  predicts 
a  nearly  correct  distance  to  G1 1.2— 0.3  (6  kpc)  and  distances 
of  9,  24,  16,  and  17,  kpc  for  Gll.4-0.1,  Gil. 15-0.71, 
Gil. 03— 0.05,  and  Gil. 18+0.11,  respectively. 

Based  on  this  “test  field,”  we  can  say  with  confidence  that 
high  surface  brightness  sensitivity,  combined  with  high  spatial 
resolution  low-frequency  radio  observations,  is  an  effective 
tool  in  the  quest  to  find  the  “missing”  Galactic  SNRs — at 
least  the  more  evolved  ones.  While  single-dish  surveys  have 
been  invaluable  for  finding  the  majority  of  SNRs  known 
today,  they  are  ineffective  at  finding  faint,  older  remnants 
because  of  confusion  with  the  Galactic  background  synchro¬ 
tron  emission,  SNRs  that  are  confused  with  thermal  sources, 
and  young  small  SNRs  for  which  they  have  insufficient 
resolution.  Interferometers  overcome  these  problems  by 
resolving  out  the  Galactic  background  while  providing  high 
surface  brightness  sensitivity  (assuming  the  inclusion  of  short 
spacing  data)  and  high  spatial  resolution.  Indeed,  with  the 
VLA  data  presented  here  we  have  more  than  doubled  the 
number  of  SNRs  in  a  ~1  deg2  field  of  view. 

However,  the  addition  of  new  SNRs  in  a  single  ~1  deg2 
field  is  insufficient  to  determine  whether  the  selection  effects 
mentioned  above  can,  in  fact,  account  for  all  the  missing 
remnants,  or  whether  there  is  a  significant  missing  link  in  our 
understanding  of  the  SNR  production  rate.  To  obtain  a  better 
statistical  sample,  we  are  currently  carrying  out  a  wider  survey 
from  /  =  +4°  to  +20°  along  the  Galactic  plane  at  330  MHz 
using  the  VLA  in  its  B,  C,  and  D  configurations.  The  resulting 
images  will  have  ^5  mJy  beam-1  sensitivity  and  ~20"  res¬ 
olution.  In  combination  with  the  planned  extension  of  the 
FIRST  survey  at  1465  MHz  to  this  region  of  the  Galactic  plane 
(R.  Becker  2003,  private  communication),  we  shoidd  be  able 
to  go  a  long  way  toward  answering  this  question. 
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